Up to eleven: radiation from particles with arbitrary energy 
falling into higher-dimensional black holes 
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We consider point particles with arbitrary energy per unit mass E that fall radially into a higher- 
dimensional, nonrotating, asymptotically flat black hole. We compute the energy and linear mo- 
mentum radiated in this process as functions of E and of the spacetime dimensionality D = n + 2 for 
n = 2, . . . , 9 (in some cases we go up to 11). We find that the total energy radiated increases with 
n for particles falling from rest (E = 1). For fixed particle energies 1 < E < 2 we show explicitly 
that the radiation has a local minimum at some critical value of n, and then it increases with n. We 
conjecture that such a minimum exists also for higher particle energies. The present point-particle 
calculation breaks down when n = 11, because then the radiated energy becomes larger than the 
particle mass. Quite interestingly, for n = 11 the radiated energy predicted by our calculation 
would also violate Hawking's area bound. This hints at a qualitative change in gravitational radi- 
ation emission for re > 11. Our results are in very good agreement with numerical simulations of 
low-energy, unequal-mass black hole collisions in D = 5 (that will be reported elsewhere) and they 
are a useful benchmark for future nonlinear evolutions of the higher-dimensional Einstein equations. 

PACS numbers: 04.70. Bw, 04.20.Dw 
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I. INTRODUCTION 

The dynamics of black holes (BHs) in generic space- 
times has attracted considerable attention in recent 
years. In astrophysics, BHs are important as sources 
of gravitational and electromagnetic waves. The inspi- 
ral and merger of BH binaries is a primary target for 
Earth-based and space-based gravitational-wave detec- 
tors In gas- rich environments, BH mergers may be 
associated with detectable electromagnetic precursors or 
afterglows @, Q and even drive the production of jets 
4\. In high-energy physics, the gauge/gravity duality 
3, [|| has created a powerful framework for the study of 
strongly coupled gauge theories, with applications in con- 
nection with the experimental program on heavy ion col- 
lisions at RHIC and LHC [H, Hj], among many others. 
The fact that BHs on the gravitational side of this corre- 
spondence are dual to thermal states of the gauge theory 
has sparked a renewed interest in BH physics. Further- 
more, some proposals to solve the hierarchy problem pos- 
tulate the existence of extra dimensions accessible only 
to gravity 8, 9]. In these scenarios, BH production from 
the collision of particles at energy scales above TeV is an 
almost inescapable consequence. 

Gravitational wave detection and high-energy applica- 
tions require an accurate knowledge of BH dynamics and 
gravitational radiation emission. This triggered research 
on the numerical evolution of the full nonlinear Einstein 
equations in four [Iol - fl8j and higher dimensions fl9l - l23j ]. 
The validation of numerical codes requires semianalytical 
tools, such as post-Newtonian theory, BH perturbation 
theory and zero-frequency expansions to model BH col- 
lisions. Such tools have been available for decades in the 
case of four-dimensional, asymptotically flat spacctimcs 



(see e.g. [24| and references therein). The same cannot 
be said of Z?-dimensional spacctimcs, but recently there 
has been significant progress in this field. For instance, 
Refs. [25|,[26[ investigated gravitational radiation and the 
quadrupole formalism in higher-dimensional, asymptot- 
ically flat spacetimes. These studies showed that odd- 
and even-dimensional spacetimes behave differently, but 
there are simple energy formulas in the Fourier-domain 
that apply to both cases [2(| . 

Linearized perturbations of higher-dimensional BHs 
arc now well understood [27]-[3(|. Historically, perturba- 
tive methods such as the close-limit approximation j32| 
(recently extended to higher dimensions [33|, HH) have 
provided guidance and insight in the numerical analysis 
of BH mergers in general relativity. The application of 
higher-dimensional BH perturbation theory to compute 
gravitational radiation in situations of physical interest 
was initiated in Ref . [3l| (henceforth Paper I) , where the 
authors studied the radiation produced by ultrarelativis- 
tic particles falling into even-dimensional, nonrotating, 
asymptotically flat BHs. 

Numerical codes to evolve the Einstein equations in 
higher dimensions are presently capable of handling low- 
energy BH collisions in five dimensions [2lj . However, 
the extension of these results to high-energy collisions in 
spacetimes of generic dimensionality presents a signifi- 
cant challenge. Motivated by these developments, here 
we extend the analysis of Paper I to study the energy and 
linear momentum radiated when particles of arbitrary en- 
ergy fall into nonrotating, higher-dimensional BHs. Our 
results are in remarkable agreement with five-dimensional 
simulations of unequal-mass BH collisions in higher di- 
mensions, that will be reported elsewhere [3^. They also 
provide useful (and sometimes surprising) insight into the 
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FIG. 1. Energy and momentum radiated plotted as functions of n for selected values of E. 



energy- and dimensionality-dependence of gravitational 
radiation produced by head-on BH collisions. 

The main findings of this paper are summarized in Fig- 
ure [TJ where we show the radiated energy (left) and lin- 
ear momentum (right) as a function of n = D — 2 for 
selected values of the particle energy per unit mass E 
(E = 1, 1.3, 1.5, 2, 3, 10). Infalls from rest correspond 
to E = 1, and the ultrarelativistic case E — >• oo is de- 
noted by "UR" in the legend. As natural in pertur- 
bation theory, the energy and angular momentum radi- 
ated are inversely proportional to the BH mass -Mbh and 
proportional to the square of the particle energy in the 
UR limit, so in the plot we normalize the radiation to 
(moE) 2 /Mbh- Paper I found that the radiated energy 
decreases with n for ultrarelativistic infalls with n < qj. 
Figure [T] shows that the total energy radiated increases 
with n for particles falling from rest. Our results for 
E = 1 and n = 3 are in remarkably good agreement 
with numerical simulations of low-energy, unequal-mass 
BH collisions in D = 5 [Hj|. They should also provide a 
useful benchmark for future nonlinear evolutions of the 
Einstein equations in higher dimensions. 

Even more interestingly, in some cases the left panel of 
Figure [T] shows the existence of a local minimum of the 
radiation as a function of n. This minimum is visible in 
the plot for the cases when the infall is not kinetic-energy 
dominated (E = 1.3 and E = 1.5), but we verified that 
it also occurs for E = 2 by extending our calculation to 
n = 11. We conjecture that such a local minimum exists 
for any E > 1, and that the radiated energy may gener- 
ically increase for sufficiently large n, eventually violat- 



ing the point-particle approximation and the area theorem 
bound. Past work showed that point-particle results in 
four dimensions can be successf ully extrapolated to the 
comparable-mass case (see e.g. [36|]). Our results imply 
that for large n this will no longer be the case. In fact, 
we find that when n = 1 1 the radiation emitted in infalls 
from rest is larger than the particle mass, and therefore 
the point-particle approximation must break down. In- 
cidentally, this breakdown fits in nicely with Hawking's 
area bound [37j . Hawking's area theorem, applied to in- 
falls from rest in generic spacetime dimension n, predicts 
that the amount of radiation emitted in equal-mass BH 
collisions decreases with n. Our results do not violate 
the area theorem. Instead, they suggest a failure of the 
point-particle approach for a change in behavior of the 
total radiated energy (in the equal-mass case) for n > 11. 
Unlike Spinal Tap's Nigel Tufncl, in higher-dimensional 
gravitational radiation we can never "go to eleven" . 

The plan of the paper is as follows. In Section |H] we 
briefly recall how to compute the radiation produced by 
particles falling radially (but with arbitrary energy) into 
a higher-dimensional BH. In Scction UHl wc present our re- 
sults on the radiated energy and linear momentum ( "kick 
velocity"), along with a preliminary comparison with nu- 
merical relativity results that will be presented in a com- 
panion paper. Section IIVI contains conclusions and pos- 
sible directions for future research. In Appendix we 
collect, for reference, some technical results. 



II. FORMULATION OF THE PROBLEM 



1 An apparent exception to this rule is the case n = 8 in Table 
VI of Paper I. Unfortunately, the extrapolated energy for n = 8 
(D = 10) was overestimated by ~ 20% in that paper. The reason 
is that we "only" computed multipoles up to I = 20 to estimate 
the total radiation, and as it turns out, this was not enough to 
get a reliable extrapolation of the total radiated energy. This 
error has been fixed here (see Table [I] below) . 



The spherically symmetric BH in D = n+2 dimensions 
is described by the Schwarzschild-Tangherlini metric (38j 
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FIG. 2. The scalar potential of Eq. (j5J as a function of the 
tortoise coordinate r, for / = 2 and selected values of n. We 
use units r h — 1. 

where dQ n is the metric of the n-dimensional unit sphere 
S n , and 
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The BH mass is related to the parameter M by 
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where A„ = 2^ n+1 ^ 2 /T[(n + l)/2] is the area of S n , 
G rl +2 is the (n + 2)-dimensional Newton constant, and c 
is the speed of light. We will set G n +2 = 1 and c = 1 in 
the following. The tortoise coordinate r* is defined by 



dr* 
dr 
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(4) 



An analytical expression for 7% (r) valid for generic n is 
given in Paper I, Eqs. (5) and (6). Here and throughout 
the paper we use the notation of Rcf . [27} ■ 

The computation of the gravitational wave emission of 
an ultrarelativistic particle plunging into a BH requires 
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the numerical integration of the inhomogencous wave 
equation for scalar gravitational perturbations ( "vector" 
and "tensor" gravitational perturbations, in the termi- 
nology of Kodama and Ishibashi, arc not excited by a 
particle in radial infall). Setting x = 2AI/r n ~ 1 , the equa- 
tion for the scalar perturbations is 
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where the scalar potential Vs is plotted in Figure [2] for 
selected values of n and 1 = 2. This potential is given by 



V s 



where the function 
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with x = 2M/r n ~ 1 , m 



f(r)Q(r) 
16r 2 H(r) 2 ' 



n(n + l)x 



1(1 



(6) 



(7) 
1), and 



Q(r) = n 4 (n + l) 2 x 3 + n(n + 1) x 

x [4(2n 2 - 'in + 4)m + n(n - 2)(n - 4)(n + l)]x 2 

- 12n[(n - 4)m + n(n + l)(n - 2)]mx 

+ 16m 3 + 4n(n + 2)m 2 . (8) 



To simplify the notation, below we will omit the su- 
perscript (n) from the wavefunction . 

Equation (O reduces to the inhomogeneous Zerilli 
equation [3^] for n = 2. The source term S 1 ;™'' in (n + 2) 
dimensions can be calculated from the stress-energy ten- 
sor of the infalling particle. Denote by E the particle 
energy per unit mass. Making use of the geodesic equa- 
tions for massive particles in radial infall 
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a straightforward generalization of the calculation pre- 
sented in Paper I yields 



S, 



00 



E 



n 2 (n + - f(r)]~2{n-2)r 
H 



V E 2 - f(r) 



(10) 



The normalized Gcgenbauer polynomials S nl are listed 
for the relevant values of n in Appendix [AJ along with 
simplified expressions of the source term in the ultrarel- 
ativistic case (E — » oo). 

We use a straightforward modification of the Fortran 
code described in Paper I to solve Eq. ((5} via Green's 
function techniques. We refer the reader to that paper 
for details. Just like in Paper I, for convenience, we set 
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the horizon radius r/j = (2M) 1 /'" -1 ) = 1 in our numeri- 
cal integrations. The energy spectrum can be expressed 
in terms of the wave amplitude at infinity given in 
Eq. (20) of Paper I, as 
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Paper I did not provide a calculation of the radiated 
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FIG. 3. Energy spectra for n = 3 and n = 6 (in units rn = !)■ 



linear momentum P*. The spectrum of the radiated mo- 
mentum can be obtained from 



dP l 
du) 



(12) 



with n l a unit radial vector on the sphere at infinity Soo . 
This results in an infinite series coupling different mul- 
tipoles. Using only the first two terms in the series, we 
find for instance 
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(13) 



(14) 



in D = 4 and D = 5, respectively. Here, <£>; denotes 
the /—pole component of the Kodama-Ishibashi wave- 
function and an asterisk denotes complex conjugation. 
We are assuming one-sided spectra. To get the total ra- 
diated linear momentum p rad ) in this work we do not 
truncate the series at the order shown in Eqs. (|13[) and 
(|14[) . Instead we sum the required number of multipolcs 
(typically ~ 10 — 20) to get the desired accuracy. 



III. RESULTS 

Our Fortran code passed several code checks. The 
spectra for n = 2 are in excellent agreement with those 
of Refs. [24|, 5(J for generic energies, and with those of 



Ref. [41( in the ultrarelativistic limit; they have been re- 
ported several times in the literature, so we do not re- 
produce them here. Our even-dimensional ultrarelativis- 
tic spectra obviously reduce to those shown in Paper I. 
Results from the Fortran code were also verified by 
comparison with a Mathematica notebook. 



A. Energy 

Figure [3] shows representative energy spectra for n = 3 
and n = 6 at different values of the particle energy. In 
the ultrarelativistic limit, as pointed out analytically in 
Ref. [25[ using Weinberg's "zero-frequency limit" approx- 
imation and confirmed numerically in Paper I, at low fre- 
quencies the spectra grow like a/ 1-2 , then they fall off ex- 
ponentially beyond a cutoff frequency u> c corresponding 
to the fundamental quasinormal mode frequency for the 
multipole in question (cf. Figure 1 in Paper I). This can 
be understood in terms of gravitational-wave scattering 
from the potential barrier surrounding the BH. The quan- 



tity oj plays the role of the energy in the Schrodinger-like 
equation (J5j) , sow 2 is equal to the maximum of the scalar 
potential Vb at first order in the WKB approximation. 
Therefore, only the radiation with energy smaller than 
the peak of the potential is backscattered to infinity; ra- 
diation with larger frequency is exponentially suppressed. 
This interpretation explains the salient features of Figure 
[31 and it is useful even in the context of comparable-mass, 
ultrarelativistic BH collisions [HJ . 

A curious new feature of the energy spectra for n > 5 
is the appearance of a double peak for large multipole 
number and intermediate particle energies. We have no 
quantitative explanation for these double peaks, but we 
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FIG. 4. Multipolar distribution for n = 3 , . . . ,9. 



suspect that they may be somehow related to the ap- 
pearance of multiple peaks in the scalar potential for low 
values of I (cf. Figure [5]). 

For a given particle energy, higher multipoles con- 
tribute more as n grows. This is even more evident when 
we look at the w-integrated multipolar components of the 
energy spectra of Figure [4] Starting from n = 6, in gen- 
eral the dominant multipole is no longer the quadrupolc. 

The total emitted energy is obtained by numerically 
integrating the spectra over w and then by summing the 
individual multipolar components AEi , which are shown 
in Figure [4] In principle, to compute the total energy we 
need to carry out a sum of all values of I up to I — > oo. 
It is of course impossible to compute multipolar contri- 
butions AEi for all values of I, so wc computed a large 
enough number of multipoles for any given dimensional- 
ity n and particle energy E. In practice, for large I we fit 
the integrated AEi with a power law of the form 

AEt = a„+ 2 r b "+ 2 , (15) 

where the coefficients (a n +2,& n +2) are obtained by fit- 
ting (typically) the last five data points of each multi- 
polar distribution in Figure U] For each n and E, the 



number of multipoles shown in the figure was chosen to 
minimize the dependence of these fits (and of the result- 
ing extrapolation) on the specific values of / chosen for 
the fit. This extrapolation introduces larger uncertain- 
ties when E and/or n get large. Our final results are 
summarized in Table U and in the left panel of Figure [1] 
For E = 1 (infall from rest) our results are well fitted 
by an expression of the form 

10 2 x ^£ rad = Cl + C2 x cf , (16) 

where a = 1.865, c 2 = 8.037x 1CT 4 and c 3 = 2.457. Now, 
based solely on the amount of emitted energy, one might 
expect the point-particle approximation to break down 
when i? rad > tjiq. Based on the extrapolation of Eq. (|16|) . 
this effectively constrains the mass ratio of the system to 
values mo/MBH < 1 when D > 13. For smaller D such 
a constraint does not apply. This may help to explain 
some results in the literature. For instance, consider the 
good agreement between numerical relativity simulations 
of equal-mass BH collisions and the point-particle extrap- 
olations to equal-mass systems. In D = 4, early work [42( 
and more recent simulations (see e.g. [HI) found that the 
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TABLE I. Total energy radiated computed using the extrap- 
olation of Eq. (|15p . For E — 2 we actually extended the 
calculation up to n = 11, and we found that a local minimum 
in the radiation occurs at n = 10: the corresponding entries 
in this table would be 0.326 (n = 10) and 0.575 (n = 11). 
For E — 1 and n = 11 the corresponding entry in this table 
would be 123, so the assumptions underlying our calculation 
are invalid (see text). 



10 2 x (M BH £ rad )/(mgff 2 ) 



n 


E = 1 


E = 1.3 


E = 1.5 


E = 2 


E = 3 


E = 10 


UR 


2 


1.04 


2.19 


3.52 


6.49 


11.9 


23.5 


26.2 


3 


1.65 


1.87 


2.75 


5.36 


9.47 


18.2 


24.9 


4 


2.02 


1.32 


1.75 


3.46 


6.48 


13.9 


19.8 


5 


2.31 


0.905 


1.00 


1.99 


4.11 


10.4 


16.5 


6 


2.92 


0.760 


0.598 


1.14 


2.67 


8.20 


13.0 


7 


4.54 


0.906 


0.457 


0.684 


1.88 


6.99 


11.4 


8 


8.27 


1.52 


0.545 


0.449 


1.44 


5.98 


10.9 


9 


17.7 


3.16 


1.00 


0.330 


1.20 


5.45 


10.6 



ratios q = mi /ma in the range between 1 and 1/4 [35j . 
Extrapolation of the numerical results to the extreme 
mass ratio limit yields [35j 

2 

E iad = 0.0164-^2_ . (is) 
Mbh 

This prediction agrees within better than 1% with the 
n = 3, E = 1 prediction listed in Table HI This excellent 
agreement provides a strong sanity check of the complex 
numerical relativity simulations, and a useful example 
of the significance of point-particle calculations such as 
those presented here. A thorough analysis of the nonlin- 
ear simulations (including more extensive comparisons 
with the point-particle limit) is in preparation [35j j . 

B. Linear momentum 



TABLE II. Momentum radiated. 



energy radiated in full nonlinear simulations of head-on 
collisions with finite mass ratio is well reproduced by sim- 
ply taking the linearized result for a particle falling into a 
BH and replacing m — > fx, where fi is the reduced mass 
of the system. This surprising agreement is at least con- 
sistent with the observation that radiation is weak, and 
therefore nonlinear effects are small. 

However, as D grows the amount of radiation emission 
also grows (a similar effect was observed recently in other 
settings [H, SI]). The extrapolation of perturbative re- 
sults to finite mass ratio must eventually break down, for 
the following reason. For two equal-mass BHs, Hawking's 
area theorem implies that the area of the final BH must 
be equal to or larger than that of the initial BHs. This 
implies the following bound on radiation emission [2~lj ] 



10 2 x (M B nP™ d )/(m 2 E 2 ) 



2M 



< 1 



BH 



(17) 



Now, this bound decreases with dimension and will even- 
tually be a strong restriction to the amount of radiation, 
violating the extrapolation of Eq. ([To]) to the equal-mass 
case. In fact, Hawking's area theorem is more restrictive 
than Eq. (|16[) for D > 13 - precisely when we know that 
the point-particle approximation breaks down due to the 
large amount of gravitational radiation. 

Summarizing, our results are self-consistent, they are 
consistent with the area theorem bound, and moreover 
they predict a qualitative change in the equal-mass BH 
collision as D increases: for D < 13 the amount of en- 
ergy lost as gravitational waves increases with D, but 
our calculation breaks down at D = 13, and presumably 
for higher dimensions the radiated energy may start to 
decrease to conform with the area bound. 

Ongoing simulations of head-on BH collisions from rest 
(E = 1) in D = 5 consider unequal-mass BHs with mass 



n 


E = 1.0 


E = 1.3 


E = 1.5 


E = 2 


E = 3 


E = 10 


UR 


2 


0.082 


0.22 


0.42 


1.1 


2.4 


5.9 


8.1 


3 


0.26 


0.25 


0.43 


1.1 


2.6 


6.8 


9.3 


4 


0.51 


0.24 


0.32 


0.82 


2.1 


6.2 


8.3 
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0.85 


0.25 


0.24 


0.59 


1.6 


5.4 


7.9 
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1.4 


0.31 


0.20 


0.43 


1.3 


4.9 


7.3 
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2.4 


0.47 


0.20 


0.31 


1.1 


4.5 


6.2 
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4.7 


0.85 


0.28 


0.23 


0.92 


4.2 


5.6 
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10 


1.8 


0.56 


0.18 


0.81 


4.1 


5.3 



In Table |TT] and in the right panel of Figure Q] we sum- 
marize the results for the linear momentum emitted in 
gravitational waves. The pattern for momentum emission 
closely mimics that of energy emission. If perturbative 
results can be extrapolated to finite mass ratios (which 
is the case for lower spacetime dimensions, see Ref. 35|) 
one expects the following mass ratio dependence [45 1: 



P = A v {m E) 



g(l-g) 
(1+g) 5 



(19) 



where q = ttiqE/Mbh- The quantity A v can be read off 
from Table ITT1 in the small-g limit. From the momentum, 
one can get the recoil velocity 



^kick _ . g 2 (l ~ g) 
C _ Av (1 + qf 



(20) 



This equation predicts a maximum kick velocity 
^kick / c = 0.0179 A„ for q = (3 + y/E)/2 ~ 0.382. Nu- 
merical simulations of BH collisions from rest in D = 5 
indicate that, in the point-particle limit [35j . 



^kick 



0.24- 



(21) 
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again in very good agreement with Table HU This is a 
nontrivial test of the simulations, because the emission of 
linear momentum involves interference between different 
multipoles. 



IV. CONCLUSIONS AND OUTLOOK 

Our results for the energy spectrum, total energy and 
momentum radiated during the head-on infall of a point 
particle into a higher-dimensional BH show an interest- 
ing and complex structure. The results indicate a beauti- 
ful concordance with the area theorem and they suggest 
that the extrapolation of perturbation theory to equal- 
mass collisions will yield wrong results for dimensions 
D > 13. This suggests that there should be a mecha- 
nism suppressing the total amount of radiation in large 
spacetime dimensions. Full nonlinear evolutions of the 
Einstein equations will probably be needed to clarify the 
exact nature of this mechanism. 

A natural and interesting generalization of our results 
would be to study the large- D limit with either numerical 
or analytical techniques. Other obvious generalizations 
include the study of infalls with finite impact parameters 
and of rotating (Myers-Perry) black holes. 

The present results should be relevant to the nascent 
field of numerical relativity in higher-dimensional space- 
times. They can be used as a guide and benchmark for 
future nonlinear simulations. Indeed, we will show in 
forthcoming work how full numerical simulations of Ein- 
stein's equations are remarkably consistent with the re- 
sults reported here f35j . 
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Appendix A: Normalization coefficients and 
ultrarelativistic limit of the source term 



for the relevant values of n and 9 = 0: 
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where X k = (I + k)(l + k- + 1). 

In the ultrarelativistic limit E — ¥ oo, the definition 
of the tortoise coordinate implies that t(r) = — r*(r). In 
units rn = (2M) 1 ^ n ~ 1 > = 1, the source term (fTU|) reduces 
to: 



(2 



S; 



s 



si 4 



s 



(5 



s. 



si 7 



s 



s 



s, 



(10 



,(11 



rm 8V4TT2 (r - l)u 2 
iwr [ U2r + 3] 2 ' 

r> 24(1 + 1) [(r 4 - r> 3 + 2(1 - r 2 )] 



(A2) 



+ 6] 2 

, r . 16 x /A 2 (2Z + 3) [(r 6 - r> 4 + 5(1 - r 3 )] 



iwr 3 V27r [i/ 4 r 3 + 10] 2 

t 20^3(2; + 4) [(r 8 - r> 5 + 9(1 - r 4 )] 



zw7rr 3 v3r [^r 4 + 15] 

r> 6VA 4 (2^ + 5) [(r 10 - r 5 )v Q + 14(1 - r 5 )] 

^rrr 4 [ Ver 5 + 21] 2 

r< 28 % /A 5 (2Z + 6) [(r 12 - r> 7 + 20(1 - r 6 )] 

V30nr [i/ 7 r 6 + 28] 2 

r> 4^A 6 (2Z + 7) [(r 14 - r> 8 + 27(1 - r 7 )] 



iumr 5 \/37r 



[v 8 r 7 + my 



18v%(2Z + 8) [(r 16 - r s )v 9 + 35(1 - r 8 )] 



i V1057 ; 



[^ 9 r 8 + 45]" 



5v/A 8 (2/ + 9) [(r 18 - r> 10 + 44(1 - r 9 )] 



Air 2 r 6 \fi 



Wior 9 + 55} z 



rm 22^(2/ + 10) 



iw37r 2 r 6 v2107rr 
[( r 20_ r i0) l/n + 54 ( 1 _ r i0)j 

[vnr 10 + 66] 2 



For the reader's convenience, here we list the normal- 
ized Gegenbauer polynomials S nl appearing in Eq. (|T0|) 



where v^. = (l + k)(l — l). These expressions are consistent 
with those listed (for even dimensions) in Paper I. 
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